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The supersymmetry method for study of disordered systems is shortly reviewed. 
The discussion starts with a historical introduction followed by an explanation of 
the idea of using Grassmann anticommuting variables for investigating disordered 
metals. After that the nonlinear supermatrix cr-model is derived. Solution of 
several problems obtained with the help of the cr-model is presented. This includes 
the problem of the level statistics in small metal grains, localization in wires and 
films, and Anderson metal-insulator transition. Calculational schemes developed 
for studying these problems form the basis of subsequent applications of the 
supersymmetry approach. 

1.1. Introduction 

The prediction of the new phenomenon of the Anderson localization^ has strongly 
stimulated both theoretical and experimental study of disordered materials. This 
work demonstrates the extraordinary intuition of the author that allowed him to 
make outstanding predictions. At the same time, one could see from that work 
that quantitative description of the disordered systems was not a simple task and 
many conclusions were based on semi-qualitative arguments. Although many inter- 
esting effects have been predicted in this way, development of theoretical methods 
for quantitative study of quantum effects in disordered systems was clearly very 
demanding. 

The most straightforward way to take into account disorder is using perturba- 
tion theory in the strength of the disorder potentialP However, the phenomenon 
of the localization is not easily seen within this method and the conventional clas- 
sical Drude formula for conductivity was considered irP as the final result for the 
dimensionality d > 1. This result is obtained after summation of diagrams with- 
out intersection of impurity lines. Diagrams with intersection of the impurity lines 
give a small contribution if the disorder potential is not strong, so that e r ^> 1, 
where £rj is the energy of the particles (Fermi energy in metals) and r in the elastic 
scattering time. 

Although there was a clear understanding that the diagrams with the inter- 
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section of the impurity lines were not small for one dimensional chains, d = 1, 
performing explicit calculations for those systems was difficult. This step has been 
done considerably later by Berezinsky^l who demonstrated localization of all states 
in ID chains by summing complicated series of the perturbation theory. This result 
confirmed the conclusion of Mott and Twose^ about the localization in such systems 
made previously. As concerns the higher dimensional systems, d > 1, the Anderson 
transition was expected at a strong disorder but it was clear that the perturbation 
theory could not be applied in that case. 

So, the classical Drude theory was considered as a justified way of the description 
of disordered metals in d > 1 and e$t 3> 1. At the same time, several results for 
disordered systems could not be understood within this simple generally accepted 
picture. 

In 1965 Gorkov and Eliashberg^ suggested a description of level statistics in 
small disordered metal particles using the random matrix theory (RMT) of Wigner- 
DysonP3 At first glance, the diagrammatic method of RefP had to work for such 
a system but one could not see any indication on how the formulae of RMT could 
be obtained diagrammatically. Of course, the description of RefP was merely a 
hypothesis and the RMT had not been used in the condensed matter before but 
nowadays it looks rather strange that this problem did not attract an attention. 

The prediction of localization in thick wires for any disorder made by Thouless^ 
could not be understood in terms of the traditional summing of the diagrams cither 
but, again, there was no attempt to clarify this disagreement. Apparently, the 
diagrammatic methods were not very widely used in that time and therefore not so 
many people were interested in resolving such problems. 

Actually, the discrepancies were not discussed in the literature until 1979, the 
year when the celebrated work by Abrahams et alP appeared. In this work, local- 
ization of all states for any disorder already in 2D was predicted. This striking result 
has attracted so much attention that it was simply unavoidable that people started 
thinking about how to confirm it diagrammatically. The only possibility could be 
that there were some diverging quantum corrections to the classical conductivity 
and soon the mechanism of such divergencies has been discovered EnHUl 

It turns out that the sum of a certain class of the diagrams with intersecting 
impurity lines diverges in the limit of small frequencies uj — > in a low dimension 
d < 2. This happens for any weak disorder and is a general phenomenon. The 
corresponding contribution is represented in Fig. 11.11 

The ladder in this diagram can be considered as an effective mode usually called 
now "cooperon" . This mode has a form of the diffusion propagator and its contri- 
bution to the conductivity a (w) can be written in the form 

a(u,)=a (l-±- [ J , ?* ) (1.1) 
y ttv J Dok 1 - iuj (2tt) J 

where Dq = Vqt/3 is the classical diffusion coefficient and ctq = 2e 2 vDo is the 
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Fig. 1.1. Diverging contribution to conductivity (coopcron) 

classical conductivity. The parameters vo and v are the Fermi velocity and density 
of states on the Fermi surface. 

Similar contributions arise also in other quantities. Eq. demonstrates that 

in the dimensions d = 0,1,2 the correction to conductivity diverges in the limit 
uj — > 0. It is very important that the dimension is determined by the geometry of the 
sample. In this sense, small disordered particles correspond to zero dimensionality, 
d = 0, and wires to d = 1. 

The contribution coming from the diffusion mode, Eq. is conceptually 

very important because it demonstrates that the traditional summation of the di- 
agrams without the intersection of the impurity lines is not necessarily applicable 
in low dimensionality. One can see that most important contributions come from 
the diffusion modes that are obtained by summation of infinite series of diagrams 
containing electron Green functions. 

The cooperon contribution, Eq. (11.11) . has a simple physical meaning. It is 
proportional to the probability for a scattered electron wave to come back and 
interfere with itselfP^ The interference implies the quantum coherence and this 
condition is achieved at low temperatures. There are many interesting effects related 
to this phenomenon but discussion of these effects and experiments is beyond the 
scope of this chapter. 

It is also relevant to mention that the cooperon contribution is cut by an external 
magnetic field, which leads to a negative magnetoresistanceP^ At the same time, 
higher order contributions can still diverge in the limit u> — > and these divergencies 
are not avoidable provided the coherence is not lost due to, e.g., inelastic processes. 

In this way, one can reconcile the hypothesis about the Wigner-Dyson level 
statistics in disordered metal particles and assertion about the localization in thick 
wires and 2D films with the perturbation theory in the disorder potential. The 
divergences due to the contribution of the diffusion modes make the perturbation 
theory inapplicable in the limit to — > and therefore one does not obtain just the 
classical conductivity using this approach. Of course, summing the divergent quan- 
tum corrections is not sufficient to prove the localization in the low dimensional 
systems and one should use additional assumptions in order to confirm the state- 
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ments. Usually, the perturbation theory is supplemented by the scaling hypothesis- 
in order to make such far going conclusions. 

At the same time, the divergence of the quantum corrections to the conductivity 
makes the direct analytical consideration very difficult for small uj because even the 
summation of all orders of the perturbation theory does not necessarily lead to the 
correct result. For example, the formulae for the level- level correlation functions^' 
contain oscillating parts that cannot be obtained in any order of the perturbation 
theory. 

All this meant that a better tool had to be invented for studying the localization 
phenomena and quantum level statistics. Analyzing the perturbation theory one 
could guess that a low energy theory explicitly describing the diffusion modes rather 
than single electrons might be an adequate method. 

The first formulation of such a theory was proposed by Wegner^ (actually, 
almost simultaneously with RefP^j). He expressed the electron Green functions in 
terms of functional integrals over conventional complex numbers S (r) , where r is the 
coordinate, and averaged over the disorder using the replica trick. Then, decoupling 
the effective interaction by an auxiliary matrix field Q he was able to integrate over 
the field S (r) and represent physical quantities of interest in terms of a functional 
integral over the N x N matrices Q, where N is the number of replicas that had to 
be put to zero at the end of the calculations. Assuming that the disorder is weak 
the integral over the eigenvalues of the matrix Q was calculated using the saddle 
point approximation. 

As a result, a field theory in a form of a so called cr-model was obtained. Working 
with this model one has to integrate over N x N matrices Q obeying the constraint 
Q 2 = 1. The cr-model is renormalizable and renormalization group equations were 
written in RefP^ These equations agreed with the perturbation theory of Eq. (jl.lj) 
and with the scaling hypothesis of RefP 

However, the saddle point approximation was not carefully worked out irP^ 
because the saddle points were in the complex plane, while the original integration 
had to be done over the real axis. This question was addressed in the subsequent 
publications E2EZI 

In the work the initial derivation of RefP^ was done more carefully shifting the 
contours of the integration into the complex plane properly. In this way, one could 
reach the saddle point and integrate over the eigenvalues of matrix Q coming to the 
constraint Q 2 — 1. After calculating this integral one is left with the integration 
over Q that can be written as 

Q = UKU-\ A =(J_° 1 ) (1-2) 

where U is an 27V x 2N pseudo-orthogonal or pseudo-unitary matrix. This matrices 
vary on a hyperboloid, which corresponds to a noncompact group of the rotations. 
This group is quite unusual for statistical physics. 
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In contrast, the method oiP^was based on representing the electron Green func- 
tions in a form of functional integrals over anticommuting Grassmann variables and 
the use of the replica trick. One could average over the disorder as well and further 
decouple the effective interaction by a gaussian integration over Q. The integration 
over the anticommuting variables leads to an integral over Q. The integral over the 
eigenvalues of Q can be calculated using, again, the saddle point method, while the 
saddle points are now on the real axis. As a result, one comes to a cr-model with 
Q-fields of the form of Eq. (|1.2|) . However, now one obtains 2N x 27V matrices U 
varying on a sphere and the group of the rotations is compact. 

The difference in the symmetry groups of the matrices Q of these two approaches 
looked rather unusual and one could only hope that in the limit N = imposed by 
the replica method the results would have agree with each other. 

This is really so for the results obtained in RefsP^^ by using the renormaliza- 
tion group method or perturbation theory. The compact replica cr-model of RefP^ 
has later been extended by FinkelsteirPDto interacting electron systems. An ad- 
ditional topological term was added to this model by PruiskerP^l for studying the 
Integer Quantum Hall Effect. So, one could hope that the replica cr-models would 
help to solve many problems in the localization theory. 

However, everything turned out to be considerably more complicated for non- 
perturbative calculations. Desperate attempts 18 to study the level-level statistics 
in a limited volume and localization in disordered wires lead the present author to 
the conclusion that the replica cr-model of RefP^ could not give any reasonable 
formulae. Calculation of the level-level correlation function using both the compact 
and noncompact replica cr-models was discussed later by Verbaarschot and Zirn- 
bauer^ with a similar result. [Recently, formulae for several correlation functions 
for the unitary ensemble (f3 = 2) have nevertheless been obtained^ from the replica 
cr-models by viewing the replica partition function as Toda Lattice and using links 
with Panleve equations.] 

The failure in performing non-perturbative calculations with the replica cr- 
models lead the present author to constructing another type of the cr-model that 
was not based on the replica trick. This method was called supersymmetry method, 
although the word "supersymmetry" is often used in field theory in a more narrow 
sense. The field theory derived for the disordered systems using this approach has 
the same form of the cr-model as the one obtained with the replica trick and all 
perturbative calculations are similar P3 

An attempt to calculate the level-level correlation function lead to a real surprise: 
the method workecP^ leading in a rather simple way to the famous formulae for 
the level-level correlation functions known in the Wigner-Dyson theoryp^ thus 
establishing the relevance of the latter to the disordered systems. Since then one 
could use the RMT for calculations of various physical quantities in mesoscopic 
systems or calculate directly using the zero-dimensional supermatrix cr-model. 

The calculation of the level correlations in small disordered systems followed by 
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the full solution of the localization problem in wiresf^l on the Bethe lattice and in 
high dimensionalityPSmn] After that it has become clear that the supersymmetry 
technique is really an efficient tool suitable for solving various problems of theory 
of disordered metals. 

By now several reviews and a book have been publisheoP^H where numerous 
problems of disordered, mesoscopic and ballistic chaotic system are considered and 
solved using the supersymmetry method. The interested reader can find all neces- 
sary references in those publications. 

The present paper is not a complete review of all the works done using the 
supersymmetry method. Instead, I describe here the main steps leading to the 
supermatrix cr-model and first problems solved using this approach. I will try to 
summarize at the end what has become clear in the last almost 30 years of the 
development and what problems await their resolution. 

1.2. Supermatrix non-linear cr-model. 

The supersymmetry method is based on using both integrals over conventional 
complex numbers Si and anticommuting Grassmann variables Xi obeying the anti- 
commutation relations 

XiXj+XjXi=0 (1-3) 

The integrals over the Grassmann variables are used following the definition given 
by BereziiP^ 



dxi = 0, / XidXi = 1 (1-4) 

With this definition one can write the Gaussian integral I a over the Grassmann 
variables as 

r N 

I A = cxp (-x + A x ) ]J dxUxi = det A, (1.5) 

i— 1 

which is different from the corresponding integral over complex numbers by presence 
of det A instead of (det .A) - in the R.H.S. In Eq. (|1.5I) . x is a vector having as 
components the anticommuting variables Xi (x + is its transpose with components 
X*) and A is an N x N matrix. 

One can introduce supervectors <£> with the components $i, 

^=( X o l ) (1-6) 



(1.7) 



A. 

and write gaussian integrals for these quantities 



N 

Is = tt~ n / exp (-$+F$) J[ d X tdXidS*dS = SDetF 
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In Eq. (jl.7[) , F is a supermatrix with block elements of the form 



Pik b. 



F ik = ( aik T ) (1.8) 



where and bik are complex numbers and Uik,Pik are Grassmann variables. The 
superdeterminant (Berezinian) SDet F in Eq. (|1.7p has the form 

SDetF = det (a- <t6~V) detfe" 1 (1.9) 

Another important operation is supertrace STr 

STrF = Tra — Trb (1.10) 

Using these definitions one can operate with supermatrices in the same way as 
with conventional matrices. Note a very important consequence of Eq. (|1.7p for 
supermatrices F that do not contain the anticommuting variables and are equal to 
unity in the superblocks Fik in Eq. f|l .8|) ( aik = bik)- In this case one obtains 

Is [Fo] = 1 (1.11) 

For such supermatrices one can write a relation that is the basis of the supersym- 
metry method in disordered metals 

Foil = J *i*ib ex P (-® + F$) d$ (1.12) 

where d$ = ir~ N fjf d X *d Xl dS*dS. 

The weight denominator in the integral in Eq. (|1.12j) is absent and this form 
is analogous to what one has using the replica trick. Applying this representation 
to correlation functions describing disordered systems one can average over the 
disorder just in the beginning before making approximations. This is what is done 
when deriving the supermatrix cr-model and let me sketch this derivation. 

Many quantities of interest can be expressed in terms of products of retarded Gf 
and advanced Gf Green functions of the Schrodinger equation. Using Eq. (I1.12j) 
one can write these functions as integrals over supervectors $ (se d 31 * 33 *) 

G?< A (y,y') = Ti [ $ a (y) (y') (1.13) 



x exp 



'/• + M(± (.-*) + «>•<*)* 



where x and y stand for both the space and spin variables. 

The Hamiltonian H in Eq. (|1.13p consists of the regular Hq and random Hi 
parts 

H = H Q + H U (H 1 )=0 (1.14) 



where the angular brackets (...) stand for the averaging over the disorder. 
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The most important contribution to such quantities as conductivity and density- 
density correlation function is expressed in terms of a product 

(r) = 2 (Gt u (r,0) G? (0, r)> (1.15) 

where r is a coordinate and w is the frequency of the external electric field. 

In order to express the function K u (r) in terms of an integral over supervectors 
one should double the size of the supervectors. Introducing such supervectors ip 
one represents the function K u (r) in terms of a gaussian integral without a weight 
denominator. This allows one to average immediately this function over the random 
part. In the case of impurities described by a white noise disorder potential u (r) 
one comes after averaging to the following expression 

K u (r) = 2 J (0) ^ (r) ^ (r ) ^ (Q) exp ^ (1 lg) 



where 



L 



dr (1.17) 



Eq. (|1.17p was obtained assuming the averages 

(u(t)u(t')) = J—6(r-r'), («(r)) = (1.18) 

where v is the density of states and r is the elastic scattering time. 

The fields tp in Eqs. (|1.16l I1.17P are conjugate to ip, the matrix A is in the space 
of the retarded-advanced Green functions and equals 

A- (J") (U.) 

The infinitesimal 6 — > +0 is added to guarantee the convergence of the integrals 
over the commuting components S of the supervectors ip. 

The Lagrangian L, Eq. (11.17[) has a form corresponding to a field theory of 
interacting particles. Of course, physically this interaction is fictitious but this 
formal analogy helps one to use approximations standard for many body theories. 

The first approximation done in the supersymmetry method is singling out 
slowly varying pairs in the interaction term. This is done writing it as 

1 f 1 

L ™t = ^ / (^) 2 dr =-j— Yl (M^) {Mp.) 

Pl+P2+P3+P4=0 



— [(V'piV'-pi+q) (^P2^- P 2-q) 



47TI/T 

Pl,P2,<?<<?0 



+ (V'piV'pJ (^-P2-q^- Pl +q) + (V'piV'ps) (V'-pi+qV'-pa-q)] l 1 ' 20 ) 

where qo is a cutoff parameter, qo < 1/1, where I is the mean free path. 

The next step is making a Hubbard-Stratonovich transformation decoupling the 
products of slowly varying pairs by auxiliary slowly varying fields. The term in the 
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second line in Eq. (|1.20[) is not important and the terms in the third line are equal 
to each other provided one uses the form of the supervectors ip of Refsl^U^ 

After the decoupling one obtains an effective Lagrangian quadratic in the fields 
ip, -0 and one can integrate out the fields ip, ip in Eq. (|1.16|) and obtain a functional 
integral over the supermatrix field Q (r) . The corresponding free energy functional 
F [Q] takes the form 



F\Q] = 



STr In 

2 



e-Hn- 



2r 



^STrQ* 
8r 



dr (1.21) 



and physical quantities should be obtained integrating correlation functions con- 
taining Q over Q with the weight exp (— F [Q]). 

The integrals with F [Q] can be simplified using the saddle point approximation. 
The position of the minimum of F [Q] is found in the limit lo — > by solving the 
equation 



Q 



i 

7T// 



H 



2t 



0(r) 



(1.22) 



One can find rather easily a coordinate independent solution of Eq. (|1.2[) . Writ- 
ing Hq in a general form as 



H = e (-iV,.) - e 



(1.23) 



and Fourier transforming the latter, one should calculate the integral over the mo- 
menta p. In the limit Eqt 3> 1 one comes to the general solution 



Q 2 = l 



(1.24) 



Although the supermatrix Q 2 is fixed by Eq. (|1.24l) . the supermatrix Q is not. 
Supermatrices Q of the form of Eq. (|1.2p are solutions for any 8x8 supermatrices 
U satisfying the condition UU — 1. With this constraint they are neither unitary 
nor pseudo-unitary as it was in RefsP^LL^ Actually, they consist of both unitary 
and pseudo-unitary sectors "glued" by the anticommuting variables. This unique 
symmetry is extremely important for basic properties of many physical quantities. 

The degeneracy of the minimum of the free energy functional F [Q] results in 
the existence of gapless in the limit 10 —> excitations (Goldstone modes). This 
are diffusion modes: so called "cooperons" and "diffusons". These modes formally 
originate from fluctuating Q obeying the constraint (|1.24[) . 

In order to write the free energy functional describing the fluctuations we assume 
that supermatrices Q (r) obeying Eq. (|1.24[) slowly vary in space. Assuming that 
to is small, lot <C 1, but finite and expanding F [Q] in this quantity and gradients 
of Q one comes to the supermatrix cr-model 



STr 



D {VQ) 2 + 2i(uj + i6)AQ dr 



(1.25) 
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where Dq = V^r/d is the classical diffusion coefficient (vq is the Fermi velocity 
and d is the dimensionality of the sample) and the 8x8 supermatrix Q obeys the 
constraint (|1.24j) . 

Calculation of, e.g., the function (r), Eq. (|1.16p . reduces to calculation of a 
functional integral over Q 

K w (r) = 2 J Q% (0) Qf a (r) exp (-F [Q]) DQ (1.26) 

Eqs. (|1.251 ll.26[) is a reformulation of the initial problem of disordered metal 
in terms of a field theory that does not contain disorder because the averaging 
over the initial disorder has already been carried out. The latter enters the theory 
through the classical diffusion coefficient Do. The supermatrix a- model, described 
by Eq. (|1.25[) resembles rj-models used for calculating contributions of spin waves 
for magnetic materials. At the same time, the noncompactness of the symmetry 
group of the supermatrices Q makes this cr-model unique. 

In order to obtain classical formulae and first quantum corrections one can 
parametrize the supermatrix Q as 



o = ir + A(i- ik 2 ) 1/2 , w=(X Q l ) {1 ' 2 '') 



and make an expansion in W in Eqs. (|1.25l ll.26j) . Keeping quadratic in W terms 
both in F [Q] and in the pre-exponential in Eq. (|1.26[) one has to compute Gaussian 
integrals over W . Fourier transforming the function K u , one obtains 

K «M = ntT ■ , (1-28) 

Eq. (jl.28l) is the classical diffusion propagator. Taking into account higher orders 
in W one can compute weak localization corrections to the diffusion coefficient. The 
first order correction is written in Eq. (II. ip . 

The precise symmetry of Q depends on the presence of magnetic or spin-orbit in- 
teractions. In analogy with symmetries of random matrix ensembles in the Wigner- 
Dyson theory^El one distinguishes between the orthogonal ensemble (both magnetic 
and spin orbit interactions are absent), unitary (magnetic interactions are present) 
and symplectic (spin-orbit interactions are present but magnetic interactions are 
absent). 

Actually more symmetry classes are possible. They are fully classified by Altland 
and ZirnbauerPsl 

In the next sections solutions of several important problems solved with the help 
of the a- model, Eq. (|1.25j) . will be presented. 

1.3. Level statistics in small metal particles. 



The first non-trivial problem solved with the supermatrix cr-model was the problem 
of describing the level statistics in small disordered metal particles. At first glance, 
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this problem is not related to the Anderson localization. However, in the language 
of the cr-model the solutions of these problems is study of the field theory, Eq. 
(|1.25p . in different dimensions. The localization can be obtained in the dimensions 
d = 1, 2 and 3, while the Wigner-Dyson level statistics can be obtained for the zero 
dimensional version of the cr-model. 

What is the zero dimensionality of the free energy functional F [Q] , Eq. (|1.25p , 
can easily be understood. In a finite volume the space harmonics are quantized. 
The lowest harmonics corresponds to the homogeneous in the space supermatrix Q. 
The energy of the first excited harmonics E\ can be estimated as 

E x = E c /A (1.29) 

where energy E c , 

E c = n 2 D a /L 2 (1.30) 

is usually called the Thouless energy. 
The other energy scale A, 

A=(vV)-\ (1.31) 

where V is the volume, is the mean level spacing. 
It is clear from Eqs. (|1.25l ll.26j) that in the limit 

E c » A, u (1.32) 

one may keep in these equations only the zero space harmonics of Q, so that this 
supermatrix does not depend on the coordinates. One can interpret this limit as 
zero dimensional one and replace the functional F [Q] by the function Fq [Q] , 

F [Q} = i ^±^STr(AQ) (1.33) 

The function R (ui) that determines the correlation between the energy levels is 
introduced as 



R(u) = ( — XI ( n ( £ k) ~ n ( £ m)) S(u- e rn + e k ) j (1.34) 

\ k.m I 

It is proportional to the probability of finding two levels at a distance to. 

Using the supersymmetry approach one can represent the functions R (w) in 
terms of a definite integral over the supermatrices Q 

R{u)= 1 -- l -Re J QllQH exp {-F [Q]) dQ (1.35) 

In order to calculate the integral in Eq. (|1.35[) one should choose a certain 
paramctrization for the supermatrix Q. 

It is convenient to write the supermatrix Q in the form 

Q = UQ U, Q =f C0S V Sin Vl> ^=(n°) (1-36) 
\ — i sin — cos 0} \ I) v I 
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where all anticommuting variables are packed in the supermatrix blocks u and v. It 
is clear that the (pseudo) unitary supermatrix U commutes with A, which drastically 
simplifies the integrand in Eq. ()1 .35|) . 

Instead of the integration in Eq. (|1.35|) over the elements of the supermatrix 
Q with the constraint (|1.24|) one can integrate over the elements of the matrix 9 
and the matrices u and v. Of course, it is necessary to write a proper Jacobian 
(Berezinian) of the transformation to these variables. The latter depends only on 
the elements of 9 and therefore the elements of u and v appear only in the pre- 
cxponential in Eq. (| 1 . 35|) . The integration over the supermatrices u and v is quite 
simple and one comes to definite integrals over the elements of 9. 

The number of the independent variables in the blocks 9 depends on the ensemble 
considered. The supermatrices Q written for the unitary ensemble have the simplest 
structure and the blocks 9 contains only 2 variables < 9 < ir and < 6\ < oo. 
The corresponding blocks 9 for the orthogonal and symplectic ensembles contain 3 
independent variables. All the transformations are described in details in RefsP^HI 

In order to get an idea about what one obtains after the integration over u and 
v in Eq. ([L36]) . I write here an expression for the unitary ensemble only 

R(u) = 1 + \ Re J exp [i (x + iS) (Ai — A)] dXidX (1.37) 

where x = ttlu/A, Ai = cosh6>i, and A = cos 9. 

So, the calculation of the level-level correlation function is reduced to an integral 
over 2 or 3 variables depending on the ensemble considered. The final result for 
the orthogonal R or th (w), unitary R un u (uj), and symplectic R sym pi (y>) ensembles 
calculated using Eq. (|1.35l) takes the following form 

_ , , sin 2 x d ( sin ir A f°° sinxt , ,„ „. 

R orth (w) = 1 — - — / — — dt (1.38) 

X (IX \ X J J -y t 

n , . sin 2 x 

R umt (w = 1 =- (1.39) 

x z 

„ , , sin 2 x d ( sin a; A f 1 sinxt , 

R sympl M = 1 - + - ) I —dt (1.40) 

Eqs. (|1.38H1.40| first obtained for the disordered metal particle^ 24 * 31 ! identically 
agree with the corresponding formulae of the Wigner-Dyson theory^' obtained from 
the ensembles of random matrices. This agreement justified the application of the 
RMT for small disordered particles suggested in RefP 

Actually, to the best of my knowledge, this was the first explicit demonstration 
that RMT could correspond to a real physical system. Its original application 
to nuclear physics was in that time phenomenological and confirmed by neither 
analytical nor numerical calculations. 

A direct derivation of Eqs. f)l -38til .4U]) from gaussian ensembles of the random 
matrices using the supermatrix approach was done in the review This allowed 



February 12, 2010 20:44 World Scientific Review Volume - 9.75in x 6.5in supersymmetry 



Anderson localization and Supersymmetry 13 

the authors to compute certain average compound-nucleus cross sections that could 
not be calculated using the standard RMT route. 

The proof of the applicability of the RMT to the disordered systems was fol- 
lowed by the conjecture of Bohigas, Giannonni and SchmicP^ about the possibil- 
ity of describing by RMT the level statistics in classically chaotic clean billiards. 
Combination of the results for clean and disordered small systems (billiards) has 
established the validity of the use of RMT in mesoscopic systems. Some researches 
use for explicit calculations methods of RMT but many others use the superma- 
trix zero-dimensional cr-model (for review see, e.g ! 34 ' 40 ' 41 l). At the same time, the 
cr-model is applicable to a broader class of systems than the Wigner-Dyson RMT 
because it can be used in higher dimensions as well. Actually, one can easily go 
beyond the zero dimensionality taking higher space harmonics in F [Q] , Eq. (|1.25l) . 
In this case, the universality of Eqs. (jl. 38111. 39)) is violated. One can study this 
limit for w> A using also the standard diagrammatic expansions of Rcf. 2 and this 
was done in RefP^l 

The other versions of the cr-model (based on the replica trick and Keldysh Green 
functions) have not shown a comparable efficiency for studying the mesoscopic sys- 
tems, although the formula for the unitary ensemble, Eq. (jl.39l) . has been obtained 
by these approaches! 20 * 43 * 

The results reviewed in this section demonstrate that the development of the 
theory of the energy level statistics in small systems and of related phenomena in 
mesoscopic systems have been tremendously influenced by the ideas of the Anderson 
localization because important results have been obtained by methods developed 
for studying the latter. 

1.4. Anderson localization in quantum wires 

The one dimensional cr-model corresponds to quantum wires. These objects are 
long samples with a finite cross-section S that should be sufficiently large, 

Sp 2 » 1, (1.41) 

where po is the Fermi momentum. In other words, the number of transversal chan- 
nels should be large. This condition allows one to neglect non-homogeneous in the 
transversal direction variations of Q. Of course, the inequality e r ^> 1 should be 
fulfilled as before. 

Then, the cr-model can be written in the form 



F[Q] 



dx, (1.42) 



where v — vS. 

Again, depending of the presence of magnetic and/or spin-orbit interactions 
the model has different symmetries (orthogonal, unitary and symplectic). It is 
important to emphasize that Eq. (|1.42p is not applicable for disordered chains or 
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thin wires where the inequality (|1.41l) is not fulfilled. However, the explicit solutions 
show that the low frequency behavior of all these systems is the same. 

Computation of the correlation function K bJ (x), Eq. (jl.26[) . with the one- 
dimensional cr-model can be performed using the transfer matrix technique. Follow- 
ing this method one reduces the calculation of the functional integral in Eq. (|1 .26|) 
to solving an effective Schrddinger equation in the space of the elements of the su- 
permatrix Q and calculating matrix elements of Q entering the pre-exponential in 
Eq. (|1.26p . This has been done in Ref.™ and presented also in the subsequent 
publications I 31 ' 33 1 

At first glance, this procedure looks very complicated due to a large number of 
the elements in the supcrmatriccs Q. Fortunately, the symmetries of the free energy 
functional F [Q] in Eq. (|1.42l) help one again to simplify the calculations. 

In order to derive the transfer matrix equations one should subdivide the wire 
into small slices and write recursive equations taking at the end the continuous limit. 
Instead of this artificial subdivision it is more instructive to consider a realistic 
model of a chain of grains coupled by tunnelling. The free energy functional Fj [Q] 
for such a chain can be written in the form 



where Jjj = J for nearest neighbors and Jij = otherwise. The summation runs 
in Eq. (|1.43|) over the grains. The coupling constant J can be expressed in terms 
of the matrix elements of the tunnelling from grain to grain Ty but at the moment 
this explicit relation is not important. 

In the limit J > 1, only small variations of the supermatrix Q in space are 
important and the functional Fj [Q] , Eq. Q1.43P , can be approximated by F [Q] , 
Eq. (|1.42[) . The classical diffusion coefficient Do corresponding to Eq. (| 1 .43[) takes 
the form 



The correlation function K u , Eq. (|1.26[) . should also be taken at the discrete 
coordinates r; numerating the grains. Then, it can be re-written identically in the 
form 



where the kernel T (ri, T2] Qi, Q2) is the partition function of the segment between 
the points r\ and ri. It is assumed that integration for this kernel is performed 
over all Q except Q\ and Q2 at the points r\ and r^. So the kernel V (r\, r%\ Qi, Q2) 
depends on supermatrices Qi, Q2 and distances r 2 — r\ (the point r 2 is to the right 
of the point ri). The function ^ (Q) is the partition function of the parts of the 




(1.43) 




(1.44) 




(1.45) 



xr (n,r a ; Qi, Qa) (Qa)™ * (Q2) dQ x dQ 2 
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wire located to the right of the point r 2 and to the left of the point rj . This function 
depends only on the supermatrix Q at the end points n or ri. 

Comparing the functions VP (Q) at neighboring grains one comes to the following 
equation 

* (Q) = In (Q, Q') Z (QO * (Q') dQ' (1.46) 



N (Q, Q') = exp (^STrQQ'J , a = 8 J (1.47) 



where 



A similar equation can be written for the kernel T (n,^; Qi, Q2) ■ Comparing this 
function at the neighboring points r and r + 1 one obtains the recurrence equation 

r (r, r'; Q, Q') - J N (Q, Q") Z (Q") T (r + 1, r'; Q", Q') dQ" 
= <5 rr ./r5 (Q - Q') (1.48) 
The (5-function entering Eq. (| 1 .48[> satisfies the usual equality 

f(Q')S(Q-Q')dQ' = f(Q) (1.49) 



Eqs. (|1.45H1.49|) reduce the problem of calculation of a functional integral over Q (r) 
to solving the integral equations and calculation of the integrals with their solutions. 
In the limit J 1 the integral equations can be reduced to differential ones. Their 
solution can be sought using again the parametrization () 1 . 36[) . The function (Q) 
is assumed to be a function of the elements of the block 8. Then, one obtains the 
differential equation for VP in the form 

•Ho* = (1.50) 

The explicit form of the operator T-Lq depends on the ensemble considered. The 
simplest equation is obtained for the unitary ensemble for which the operator Hq 
takes the form 



Ho — — : 



2-kvDq 
where 



1 9 J 9 + 1 9 J 9 



Ja = (Ai - A)" 2 



i (ui + iS) ttv (Ai — A) (1.51) 



Similar equations can be written for the central part entering Eq. (|1.45p . 

Solving these equations and substituting the solutions into Eq. (|1 .45|) one can de- 
termine (at least numerically) the frequency dependence of the function (n, r^) 
and, hence, of the conductivity for all frequencies in the region ujt <C 1 and distances 
In. - ri\po » 1. 
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The calculation becomes considerably simpler in the most interesting case of low 
frequencies ui <C {y 2 Dq) . In this limit, the main contribution into the integral in 
Eq. (jl.45j) comes from large Ai 3> 1 and the solution \I> of Eq. (|1.50j) is a function 
of only this variable. 

Introducing a new variable 

z = -icj27r 2 £ 2 L>oAi (1.52) 
one can reduce Eq. (|1 .51|) to the form 

d 2 ^ (z) 

-z^P + *(z) = (1.53) 

with the boundary condition 

*(0) = 1 (1.54) 
The Fourier transformed function K u (k) takes the form 

K u (k) = r^, A(k)= (<S> k (z) + <f>_ k (z))<f(z)dz, (1.55) 



where the function & k (z) satisfies the following equation 
d ( 2 d<S> k (z) 



,'_2_w ) +ikL c <S> k (z) + z$ k (z) = V{z) (1.56) 
dz \ dz J 

with the length L c equal to 

L c = 2-kvSD (1.57) 

The length L c is actually the localization length, which will be seen from the final 
result. Equations (|1 .53111 .5"6")) can also be obtained for the orthogonal and symplectic 
ensembles but with different localization lengths L c . The result can be written as 

j^symplectic ^j^unitary ^j^orthogonal /-^ 

The residue of the function K u is proportional to the function (r, r', e) introduced 
by AndersonfS 

Poo (r, r', s) = J2 \<h Ml' \<t>k (r')\ 2 6 (e - e k ) , (1.59) 

k 

where (j) k (r) are exact eigenfunctions. 

Eqs. (|1.53H1.5B|l exactly coincide with the low frequency limit of equations de- 
rived by Berezinsky^ provided the length L c is replaced by the mean free path I, 
which shows that the low frequency limit of the one dimensional systems is universal. 

The exact solution of Eqs. (|1.53II1.5B| leads to the following expression 
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In the limit x 3> L c , Eq. (|1.60p reduces to a simpler form 



The exponential form of p^ (x) proves the localization of the wave functions and 
shows that the length L c is the localization length. Note, however, the presence of 
the pre-exponential |x| 3 ^ 2 . Due to the factor the integral over x of Poo (x) remains 
finite even in the limit L c — > oo. Actually, one obtains 



Poo (x) dx = v, (1-62) 

which proves the localization of all states. 

At small k <C L^ 1 , the function A (k) in Eq. (| 1 . 55[) takes the form 

A(k) = l-4((3)k 2 L 2 c 



and the static dielectric permeability e equals 
e = —Aixe 2 v 



, d 2 A(k) 



dk 2 



= 32C(3)e A vL A c (1.63) 

fc=0 



where £ (x) is the Riemann C-function. 

All these calculations have been performed for a finite frequency u> and the 
infinite length of the sample. One can also consider the case of the zero frequency 
and a finite length L. A full analysis of this limit has been presented by Zirnbauei^l 
who calculated the average conductivity as a function of L. 

There is another Fokker-Planck approach to study transport of disordered wires 
developed by Dorokhov, Mello, Pereyra, and Kumar^^ES (DMPK method). It can 
be applied also to thin wires with a small number of channels. At the same time, this 
method cannot be used for finite frequencies. In the case of thick wires with a large 
number of the channels and zero frequencies, the equivalence of the supersymmetry 
to the DMPK method has been demonstrated by Brouwer and FrahmPsl 

Many interesting problems of banded random matrices^! and quantum chaos 
(like kicked rotor^) can be mapped onto the ID supermatrix cr-model. However, 
a detailed review of these interesting directions of research is beyond the scope of 
this paper. 

1.5. Anderson localization in 2 and 2 + e dimensions. 

Study of localization in 2 and 2 + e using the replica cr-model was started by Weg- 
nei^ using a renormalization group (RG) technique. He was able to write the RG 
equations for the orthogonal and unitary ensembles that could be used in 2 dimen- 
sions and extended into 2 + e dimensions for e <C 1. The latter was done with a 
hope that putting e = 1 at the end of the calculations one could extract at least 
qualitatively an information about the Anderson metal-insulator transition in 3 di- 
mensions. Based on this calculation a conclusion about the localization at any weak 
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disorder in 2D was made. As concerns 2 + e, an unstable fixed point was found, 
which following the standard arguments by Polyakov- 5 - 1 signaled the existence of the 
metal-insulator transition. 

The symplectic case was considered within the compact replica cr-model in RefP^ 
using the same method of RG and it was shown that the resistivity had to vanish in 
the limit of ui — > 0. The difference between the replica cr-models used in Refs LUCIUS 
(noncompact) and RefP^ (compact) is not essential when applying the RG scheme. 

Exactly the same results are obtained with the supermatrix cr-model using the 
RG techniqu e 1 31 * 33 * 52 ! and let me sketch the derivation here. As usual in the RG 
method, one introduces a running cutoff parameter and coupling constants depend- 
ing on this cutoff. The cr-model for such couplings can be written as 

F = j J STr (VQ) 2 + 2iuAQ dr (1.64) 

where Co = uj /Do- The bare value of t equals t = 8 (ttj/D ) (c.f. Eq. (pL25j> ). 

The cr-model looks similar to classical spin cr-models considered in Refp^J and 
one can follow the RG procedure suggested in that work. Using the constraint 
(|1.24j) one can write the supermatrix Q in the form 

Q = VAV, (1.65) 

where VV = 1 so that V is a pseudo-unitary supermatrix. 

In order to integrate over a momentum shell one can represent the supermatrix 
V in the form 

V{r)=V (r) V (r) , (1.66) 

where Vq is a supermatrix fast varying in space and V is slowly varying one. These 
supermatrices have the same symmetry as the supermatrix V. 

Substituting Eq. (|1.66[) into Eq. (|1.64p one can write the free energy functional 
F \Q] in the form 



F = 



l - J STr [(VQ ) 2 + 2 [Q , VQ ] $ + [Qo, $] 2 + 2iuVAVQ ] dr (1.67) 



Q = VqAVq, $ = VWV = 

The next step of the RG procedure is to integrate over the fast varying matrices 
Qo and reduce to a functional containing only slowly varying variables V. After 
this integration the free energy F in Eq. (|1.67|) should be replaced by energy F 
describing the slow fluctuations 

F = — In J exp (-F) DQ (1.68) 

The integration over the supermatrix Qo can be done using a parametrization (|1.27[) 
or a more convenient parametrization 

Qo = A(l + P)(l -Py\ PA + AP = Q. (1.69) 
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Integration over the fast variation means that one integrates over Fourier trans- 
formed Pk with A/co < k < ko, where ko is the upper cutoff and A < 1. As a result 
of the integration one comes to the same form of the functional F as in Eq. (| 1 .64[) . 
The constant uj does not change under the renormalization but the new coupling 
constant t can be written as 

V 8A fco fc 2 (2^)V 1 ' 

The correction to the coupling constant t, Eq. (|1.70p . is written in the first order 
in t. The parameter a depends on the ensemble and equals 

!— 1, orthogonal 
0, unitary (1-71) 
1 symplectic 

Stretching the coordinates in the standard way and changing the notation for the 
coupling constant t — > 2 d+1 ndT (d/2) t, where T is the Euler T-function one obtains 
the RG equation for t 

ft(t) = ^ = (d-2)t + at 2 (1.72) 

where ft (i) means the Gell-Mann-Low function. 

In 2D, the solution of this equation for the coupling constant t (proportional to 
resistivity) takes the form 

t (u) = r , . . (1.73) 

For sufficiently high frequencies uj the resistivity and the diffusion coefficient D (uj) 
proportional to i _1 (uj) coincide with their bare values. 

Decreasing the frequency ui results in growing the resistivity for the orthogonal 
ensemble until the coupling constant t (uj) becomes of the order 1. Then, the RG 
scheme is no longer valid because the expansion in t in the R.H.S. of Eq. (|1.72[) is 
applicable only for t -C 1. However, it is generally believed that t diverges in the 
limit uj — > and this should mean the localization of all states with an exponentially 
large localization length 

L c cxexp(l/i ) (1.74) 

In the symplectic ensemble the resistivity t (w) decreases with decreasing the 
frequency u>. This interesting result was obtained in the first order in to by Hikami, 
Larkin and Nagaoka. 53 However, Eq. (I1.73|) means moreP^ If the bare to is small, 
to <C 1, the effective resistivity t (uj) decays down to zero in the limit u> — > 0. In this 
case the constant t (uj) is small for any frequency and the one loop approximation 
used in the derivation of Eq. (|1.72|) is valid for all frequencies. So, the solution for 
the symplectic ensembles, when used for the low frequencies, is the most reliable 
one obtained with the RG method. 
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As concerns the unitary ensemble, the first order contribution vanishes and one 
should calculate corrections of the second order. As a result, one comes to the 
following dependence of t (w) on the frequency 

t H = - m (1.75) 

One can see from Eq. (11.75)) that the resistivity t (lj) grows, as in the orthogonal 
ensemble, until it becomes of order 1. Again, this behavior is interpreted as lo- 
calization for any disorder. The conclusions about the localization in 2D for the 
orthogonal and unitary ensembles were made first in RefP^ and this agreed with 
the results based on using the scaling hypothesisP 

Wegner developed also theory of the Anderson metal-insulator transition in the 
dimensionality 2 + e for e < lN' One can see that the RG equation (|1.72[) has a 
fixed point t c — e, at which the Gell-Mann-Low function vanishes. At this point the 
total resistance of the sample does not depend on the sample size and this point 
should correspond to the Anderson metal-insulator transition. 

Linearizing function (3 (t) near the fixed point t c one can solve Eq. (|1.72[) . As a 
result one can find a characteristic (correlation) length £ near the fixed point 

^^(^ir) 1/V ' y= -p ,{tc) (L76) 

where £o is the size of a sample having the entire resistance to. Assuming that the 
length £ is the only characteristic length in the system and that the conductivity a 
is proportional to t~ 1 £ 2 ~ d , one can write the equation for the conductivity in the 
following form 

" 2 ^c-t V _ d-2 

The explicit values of the critical resistance t c and the exponent s for the or- 
thogonal and unitary ensembles equals 



- = A-^-[^) , s=— (1.77) 



d — 2, orthogonal 
(2(d-2)) 1/2 , unitary (1 ' ?8) 



and 

1 



1/2 ( L7Q ) 

demonstrate that the metal-insulator transition exists in any 
dimensionality d > 2 and the conductivity near the transition obeys a power law. 
Of course, this consideration is restricted by small e = d — 2 and one can use the 
result in 3D only qualitatively. 

The scaling approach developed for small e is similar to the one developed for 
conventional phase transitions in, e.g., spin models where one can also write a- 
models. This method is not sensitive to whether the symmetry of the supermatrices 
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Q is compact or noncompact. Using this approach one comes to the conclusion that 
the Anderson metal-insulator transition is very similar to standard second order 
phase transitions. 

In the next section the same problem will be considered on the Bethe lattice 
or in a high dimensionality. Surprisingly, the result will be very different and the 
peculiarity of the solution originates from the noncompactness of the group of the 
symmetry of the supermatrices Q. 

1.6. Anderson metal-insulator transition on the Bethe lattice or in 
a high dimensionality. 

It is generally difficult to find the critical point for a transition between differ- 
ent states and describe the critical behavior in its vicinity. The Anderson metal- 
insulator transition is definitely not an exception in this respect. Usually, identifying 
a proper order parameter one can get an idea about a transition using a mean field 
approximation. As concerns the Anderson transition, this is not possible. Although 
the a- model, Eq. (|1.25|) . looks very similar to spin models in a magnetic field, one 
cannot take an average of Q with the free energy F [Q] , Eq. (|1.25j) . as the order 
parameter because it determines the average density of states and is not related to 
the Anderson transition. 

At the same time, the mean field approximation works very well in high dimen- 
sionality or on special structures like the Bethe lattice. 

The Anderson model of the Bethe lattice was studied for the first time by Abou- 
Chacra, Anderson and Thoulessf^ who proved the existence of the metal-insulator 
transition and found the position of the mobility edge. With the development of 
the supersymmetry technique it became possible to describe the critical behavior 
both in the metallic and insulating regime. Considering a granular model one could 
obtain results for the orthogonal, unitary and symplectic ensembles. Later the 
Anderson model has also been described. 

It turned out that in all the cases the critical behavior was the same, which 
contrasts the results obtained within the 2 + e expansion. This could not be a big 
surprise because for most phase transitions the high dimensional results are more 
"universal" than those obtained in lower dimensions. However, the results for the 
metallic and insulating regimes did not obey the conventional scaling and this was 
completely unexpected. 

The first attempt to solve the granular version of the supermatrix a-model on 
the Bethe lattice has been undertaken in RefpS In this work correct integral equa- 
tion have been written for description of critical behavior near the metal-insulator 
transition and the position of the mobility edge has been found. However, attempts 
to find a solution of this equation related to scaling properties of the 2 + e limit 
were not successful, which lead to wrong conclusions. 

Studying numerically the integral equation derived irP^l Zirnbauer^S found a very 



February 12, 2010 20:44 World Scientific Review Volume - 9.75in x 6.5in supersymmetry 



22 K.B. Efetov 

unusual behavior near the critical point and presented formal reasons explaining 
this behavior. Finally, the density-density correlation function has been calculated 
for the unitary^ and orthogonal and symplectic ensembles.^ This determined 
the diffusion coefficient in the metallic region and localization length and dielectric 
permeability in the insulating one. 

The form of the density-density correlation function on the Bethe lattice dif- 
fers from the one on conventional lattices. Therefore the problem of the Anderson 
localization has been considered on such lattices in an effective medium approxima- 
tion.^ The latter becomes exact on the real lattices in a high dimensionality d^> 1 
and the basic equations and results are similar. The derivation of the equations and 
the final results are shortly displayed below. A detailed discussion can be found 
irP 

The scheme of the derivation of the equations is similar to the one presented 
in Sec. 11.41 for one dimensional structures consisting of the grains. We start with 
Eq. (| 1 .43[) written on a <i-dimensional lattice with d 3> 1 or on the Bethe lattice. 
Denoting by ^ (Q) the partition function of a branch of the tree structure with 
a fixed value Q at the base and comparing it with the partition function on the 
neighboring site one comes to a non-linear integral equation 

*(Q) = J N (Q, Q') Z a (QO * m (QO dQ' (1.80) 

where m — 2d — 1 for a (i-dimensional lattice and is the branching number on the 
Bethe lattice. The functions N (Q, Q') and Zq (Q) have been introduced in Eq. 

The case m = 1 corresponds to the one-dimensional chains of the grains and Eq. 
(|1.80p coincides with Eq. (| 1 .46[) in this limit. In this particular case equation (|1.80p 
is linear and, as we have seen in Sec. 11.41 all states are localized for any disorder. 
However, at m > 1 the integral equation (jl.80l) is non-linear and has a bifurcation 
at a critical a c corresponding to the Anderson metal-insulator transition. 

The density-density correlation function K u , Eq. (| 1 . 26|) . can be written in the 
form 

K u (n, r 2 ) = -2vrW J Q^P 33 (r, Q) Z (Q) * (Q) dQ (1.81) 

where the function P (r, Q) satisfies for the high dimensional lattices the following 
equation 

P(r,Q)-J2w(r- r') f N (Q, Q') P (r', Q') Z (Q') dQ' (1.82) 

+ m j ' N 2 (Q, Q') P (r, Q') Z (Q 1 ) dQ' = S (r) Q 21 * (Q) . 
In Eq. (|1.44|) the function N 2 (Q, Q') is equal to 

^2 (Q, Q')= In (Q, Q") N (Q", Q) Z (Q") dQ" 
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and 

W(r-r') = {^ |r - r ' l = 1 
1 ' \0,|r-r'|^l 

The third term in the L.H.S. of Eq. (|1.82j) takes into account the fact that two 
segments of a broken line cannot coincide. Eqs. (jl.81l ll.82j) are very similar to 
Eqs. ([1.451 11.481) written for the ID case. This is natural because in both the cases 
loops are absent. Their solution for a function ^> (Q) found from Eq. (|1.80j) can be 
obtained making a spectral expansion of P (r, Q) in eigenfunctions of the integral 
operators entering the L.H.S. of Eq. (|1.82[) . 

In principle, this procedure is straightforward. However, solving the integral 
equation (|1.80|) is not simple because it contains a large number of the elements of 
the supermatrix Q. 

Fortunately, Eqs. (|1 .80111 .82)) drastically simplify in the metallic regime near the 
metal-insulator transition and everywhere in the insulating regime provided one 
considers the low frequency limit u> — > 0. The formal reason for this simplification 
is that the main contribution into the correlation functions comes in these cases 
from the region of very large values of the variables Ai > A/cj ^ 1. The same 
simplification has helped one to solve the problem of the localization in wires in 
Sec. O 

Nevertheless, the full analysis is quite involved even for small oj. Details can be 
found again in Ref.ESJ and here I display only the final results. 

In the insulating regime, a < a c , only = 1 is the solution of Eq. (| 1 .80[) in the 
limit oj = 0. This solution of the simplified equation persists for all a but another 
solution appears in the region a > a c . The latter solution considered as a function 
of 9\ = ln(2Ai) has a form of a kink moving to infinity as a — » a c . The position 
9\ c of the kink depends on the distance from the critical point a c as 

9 lc = s(a-a c y 1/2 (1.83) 



where s is a number of order 1 . The dependence of \I/ (Ai ) is represented in Fig. 11.21 
Only this solution should be used for a > a c and this leads to a very non-trivial 
critical behavior of the diffusion coefficient. 

The position of the critical point a c and the critical behavior have been calcu- 
lated for all 3 ensembles. For large m, the value a c for the orthogonal and unitary 
ensembles is determined by the following equations 

(§£ ) V2 m ln £ = !. orthogonal 
V2ttJ TOm ~ = ' unitary 

One can see from Eq. (|1.84|) that the metallic region is broader for systems with 
the broken time reversal invariance. In other words, applying a magnetic field shifts 
the metal- insulating transition to larger values of a c . This result correlates with the 
one, Eq. ()1.78[) . obtained in 2 + e dimensions. 
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Fig. 1.2. Numerical solution ^(Ai) for m = 2 and different hopping amplitudes in the critical 
metallic regime. The inset shows 9jj defined by ^(coshdjj) = 0.5 as a function of (a — etc) -1 / 2 

Although the position of the Anderson transition depends on the ensemble con- 
sidered, the form of the correlation functions is the same. 

In the insulating regime, the function poo (r), Eq. (|1.59[) . takes for r ^> L c the 
following form 

Poo ( r ) = const r-^+V^L-^exp f~^-) ( L85 ) 

where L c is the localization length. 

Near the transition the localization length L c grows in a power law 

const 



L r . = 



(1.86) 



[a c - a) 

In this regime there is another interesting region of 1 -C r <C L c where the function 
Poo (t) decays in a power law 

Poo ( r ) = const r^ 11 - 1 (1.87) 



Remarkably, Eq. (|1.85l) obtained for d ^> 1 properly describes also the one- 
dimensional wires (c.f. Eq. (ILGip - ). 

The integral of Poo (r) over the volume is convergent for all a < a c and remains 
finite in the limit a — > a c indicating that the wave functions at the transition point 
decay rather fast. At the same time, all moments of this quantity diverge in this 
limit. The second moment determines the electric susceptibility re, 

nS a/ 3 = e 2 J Tarppaz (r) d d r (1.88) 

Near the transition calculation of the integral in Eq. (|1.88p leads to the result 

k = 4it 2 vcL c (1.89) 

where c is a coefficient. 
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This equation shows that the susceptibility in the critical region is proportional 
to the localization length L c and not to L\ as it would follow from the one-parameter 
scaling^ and results obtained in 2 + e dimensions. The unusual dependence of the 
susceptibility k on L c in Eq. (|1.89[) arises formally from the anomalous exponent 
(d + 2) /2 in the power law behavior of the pre-exponent in Eq. (jl.85j) . 

As concerns the metallic regime one comes for the real lattices to the diffusion 
propagator, Eq. (|1.28[) . at all a > a c . However, except the limit a ^> 1, the diffusion 
coefficient D obtained now is different form the classical diffusion coefficient Dq. Its 
behavior in the critical region a — a c <C a c is especially interesting. This is not a 
power law behavior as one could expect from the one parameter scaling. Instead, 
the diffusion coefficient decays near the transition exponentially 



exp 

D = const 



I \-l/2' 

-s(a-a c ) 



(1.90) 



i ^3/2 
(a - a c ) 

This is a very unusual behavior. Formally, it follows from the non-compact sym- 
metry of the supermatrices Q. For any compact symmetry one would obtain in the 
same approximation a power law dependence of the diffusion coefficient on a — a c . 
The exponential decay of the diffusion coefficient D, Eq. (|1 .90|) . follows from the 
shape of the function 'J, Fig. 11.21 The position of the kink Xh = cosh^, Eq. 
(|1.86|) . goes to infinity as Xh oc exp s {a ~ Uc) 1 ^ 2 and this results in the form 
(|1.90[) of the diffusion coefficient. 

The same results, Eqs. (II .85111 .901) . have been obtained latei^for the Anderson 
model on the Bcthc lattice and this completed the study of this model started in 
Ref.^ The agreement of the results obtained for the Anderson model and granulated 
cr-model is, of course, not accidental because, the critical behavior is formed by long 
time correlations and the result should not be sensitive to short distance structures. 
As it has been discussed previously, the low frequency behavior of wires and strictly 
one dimensional chains is also described by identical equations. 

The exponential decay of the diffusion coefficient was interpreted in RefP^ in 
terms of tunnelling between quasi-localized states. This may happen provided the 
wave function is concentrated in centers with a large distance 

C oc (a- a c y 1/2 (1.91) 

between them. The decay of the amplitudes of the wave functions in the single 
center is fast as it can be seen from the fast decay in Eq. (|1.87[) . Then, the 
tunnelling leads to an overlap between the wave functions of the different centers 
and to formation of a conduction band with an effective bandwidth T, 

rcxAexp(-aC) (1.92) 

where a is a coefficient. 

The exponential decay of the diffusion coefficient D, Eq. (|1.90j) . can follow 
quite naturally from such a picture. Of course, the picture implies the existence of 
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weakly overlapping centers of the localization near this transition. A well established 
multifractality of wave functions at the transition (for a recent review, see, e.gP^) 
may point out on this strong inhomogeneity near the transition. 

Another indication in favor of the presented picture comes from the fact that 
the solution ^> of Eq. (jl.80j) looses the sensitivity to the existence of the transition 
at frequencies u > T. This can be seen from a more detailed analysis of Eq. (|1.80p . 
The interpretation in terms of formation of a very narrow conduction band near the 
transition with the bandwidth V is consistent with this property of the solution it . 

The fixed point found in 2 + e for small e corresponds to a weak disorder and 
the strong inhomogeneities are not seen in this approach. One cannot speak of a 
narrow conduction band near the transition in 2 + e dimensions within this picture. 

In principle, centers of (quasi) localization exist in 2D and can be described in 
the framework of er-modcl (for a review, see RefsP^^J). However, the idea about 
these centers of the (quasi) localizations is not incorporated in the conventional 2 + e 
scheme. So, the standard continuation of the results obtained for small e to e = 1 
may result in loosing an important information. 

The non-trivial form of the function ^ (see e.g. Fig. I1.2[) has lead the present 
author to the idea22l that this function might play the role of an order parameter for 
the Anderson transition. It was guessed that a Laplace transform of this function 
could be related to a conductance distribution. This idea has been further developed 
in Ref.f^ where a functional in an extended space was constructed such that its 
minimum was reached at the function ^ (Q). This resembles the Landau theory of 
phase transitions but the role of the order parameter is played a by a function. 

The concept of the function order parameter was also discussed in later works 
on the Bethe latticed 

1.7. Discussion. 

In this paper the basics of the supersymmetry method has been presented. It 
is explained how the non-linear supermatrix er-model is derived and it is shown 
how one can calculate within this model. It is demonstrated how one comes to 
the Wigner-Dyson statistics in a limited volume and how one obtains Anderson 
localization in disordered wires. Renormalization group scheme is explained in 2 
and 2 + e dimensions for small e, renormalization group equations are written and 
solved. It is shown how one solves the problem of the Anderson metal-insulator 
transition on the Bethe lattice and high dimensionality. 

From the technical point of view all this was a demonstration how one can 
calculate in the dimensions d = 0, d = 1, d = 2, and d > 1. Due to the lack 
of the space the present paper is not a complete review of the application of the 
supersymmetry technique and many interesting works are not mentioned. However, 
the calculational schemes presented here have been used in most of the subsequent 
works. So, having read this paper one can get an idea on how one can work in all 
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situations where the supersymmetry method is useful. 

This is a chapter in the book devoted to 50 years of the Anderson localization 
and I tried to describe shortly how one of the directions of the field was developing in 
1980s after the second most important work on the Anderson localization^ has been 
published. Many of the authors of the present volume entered this field motivated 
by this publication. I hope that the development of the supersymmetry method has 
been useful in solving several interesting problems of the Anderson localization. 

Although the supersymmetry method proved to be an adequate method for 
studying disordered systems (at least, without electron-electron interaction), several 
very important problems have not been solved so far. In spite of the common 
believe that all states are localized in disordered films (orthogonal and unitary 
ensembles) , the solution for the two-dimensional cr-model has not been found in the 
limit of low frequencies. The problem of the integer quantum Hall effect has not 
been solved either, although the idea about instantons^l was very useful for the 
understanding of this phenomenon. The problem of describing the critical behavior 
near the transition between the Hall plateaus still awaits its resolution. 

One more interesting problem is to understand the critical behavior near the 
Anderson transition. 

Of course, a lot of information comes from numerical simulations but solving 
the 2D problem analytically would be really a great achievement. As concerns the 
Anderson transition in 3D, the hope to solve it exactly is not realistic because even 
a simpler Ising model has not been solved in spite of numerous attempts. However, 
in the conventional theory of phase transition one can start with a mean field theory 
justifiable in high dimensions, determine the upper critical dimension and then make 
an expansion near this dimensionality. 

Unfortunately, until now a similar procedure has not been found for the Ander- 
son transition, although the supermatrix cr-model resembles spin models for which 
this procedure is standard. This concerns also the 2D case, where conventional spin 
cr-models are solvable. However, the well developed methods like the Bethe Ansatz 
or methods of the conformal field theory do not work here. 

The formal reason of the failure of these approaches for studying the supermatrix 
cr-model is that the group of the symmetry of the supermatrices Q is not compact. 
These supermatrices consist of a block varying on a sphere and another one with 
elements on the hyperboloid. The latter part of Q is formally responsible for the 
localization but its presence leads to difficulties when applying the well developed 
methods. It is clear that the importance of the noncompact symmetry is not fully 
appreciated. 

I can only express my hope that these problems will be resolved in the next 50 
years and the book devoted to 100 years of the Anderson localization will contain 
the complete theory of this phenomenon. 
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